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Abstract. Probably the most fundamental characteristic of the quasar- 
AGN power house, the mass of the central black hole, is the least well 
known. The broad emission lines have been, and probably will remain, 
our best probe of the central mass. Using these probes to estimate 
the black hole mass suggests that over more than six orders of magni¬ 
tude, the ratio between the continuum luminosity and the central mass 
(the Eddington ratio) has a spread perhaps as narrow as 1-2 orders 
of magnitude, while other methods give a larger spread and possibly a 
luminosity-dependent Eddington ratio. I review the three main classes 
of mass estimation methods—BLR kinematics. X-ray variability and ac¬ 
cretion disk modeling, and their results for L/LEdd- Potential sources of 
error and biases are discussed. 


1. The L/M relation and mass estimation methods 

The most fundamental, yet the least well known, parameter of the AGN central 
engine is its mass. While the AGN luminosity, which is readily measurable, 
varies over six orders of magnitude, the ratio L/M (or, alternatively, the Ed¬ 
dington ratio, L/Le) apparently has a rather small range, indicating that the 
L/M value is a property characterizing the AGN phenomenon. In terms of the 
Eddington ratio we find from various mass estimation methods and for various 
AGN samples L/Le ~ 0.001-1. This value and its spread depend on the method 
used to estimate the mass. 

The methods used to determine the central mass in AGN can be divided into 
three major groups: BLR kinematics, X-ray variability and accretion disk spec¬ 
tral fitting. The first group uses the broad emission lines as probes, assuming 
the velocity dispersion of the line-emitting gas is induced by the gravitational 
potential of the central mass. The second gives the variability mass limit (eq. 

, while the third is a variant of the temperature mass (eq. , essentially trying 
to fit the UV bump with a spectrum from a thin accretion disk model. In the 
following sections I will describe the three groups and the problems associated 
with each of them. 


2. Some fiducial masses 

Any compact, accretion-powered radiation source can be assigned several fiducial 
mass estimates: 
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1. The Eddington limit. In order to maintain steady spherical accre¬ 
tion the luminosity must be less than the Eddington luminosity, L < L^dd = 
AnGMuipc/aT = 1-3 x 10^®Mg erg s“^, or Mg > L 46 where Mg = M/IO^Mq 
and L 46 = L/ 10 ^® erg s“^. 

2. The variability limit. The shortest time scale for global variations in 
the luminosity is the light travel time across the Schwarzschild radius, Rs/c = 
2GM/c^ = 10^Mg s. Hence if the luminosity is observed to vary significantly on 
a time scale At, the black hole mass has to be 

Mg < (At/lO^-^ s), (1) 

where r is the effective radius of emission in units of Rg- 

3. The temperature mass. If a luminosity L comes from an accretion disk re¬ 
gion of radius R and temperature T, then L < AnR'^aT'^. If a spectral feature 
at frequency u is due to emission from a black body accretion disk at a charac¬ 
teristic radius R, then the temperature is given by hi^ = 3kT, and we have an 
upper limit on the black hole mass; 

Mg > 100(r/105 K)-^Lli\R/Rg)-\ (2) 

4. The evolutionary mass. This mass involves the evolution of a massive 
black hole due to accretion. We may define the Eddington time as the time 
required for an accretion-powered object radiating at the Eddington luminosity 
to double its mass, tE = M(P‘ /Le = 4 x 10 ®y. The observed luminosity implies 
an accretion rate of M = O.16e“^L46M0y“^ where e is the efficiency. Combining 
these two expressions gives the evolutionary mass, accumulated if the accretion 
rate is maintained during a time tE ■ 

Mg = 0.6e“^r?“VL-^46 (3) 

where the parameters rj and fi are the Eddington luminosity ratio and the 
fraction of the active time for intermittently active AGN. 

3. BLR kinematics 

The broad emission lines characteristic of quasar and Seyfert 1 spectra are be¬ 
lieved to originate from gas highly ionized by the continuum radiation of the 
central source. The line width must be due to bulk motions, probably an en¬ 
semble of clouds moving in the gravitational potential of the central black hole . 
Alternative scenarios include radiative acceleration (e.g. Blumenthal and Math¬ 
ews 1975), which would induce radial outflow. However, radial motions are 
probably excluded because they would induce asymmetric line profiles (if the 
emitting clouds are optically thick) which are usually not observed. Most work¬ 
ers today agree that the line width in AGN is induced by Keplerian bulk motions. 
If this is the case, knowledge of the distance of the line-emitting gas from the 
central mass would allow us to estimate the black hole mass, 

M v^R/G. (4) 

The kinematic mass-estimation methods can be divided into two sub-groups, de¬ 
pending on the method used to find the emission-line distance: photoionization 
methods and reverberation mapping. 
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3.1. Photoionization: the ionization-parameter method 


We assume that the spatial extent of the BLR may be represented by a charac¬ 
teristic size — e.g. the radius at which the emission peaks. Indeed, reverberation 
analyses have shown that the time lag of emission-line variability varies depend¬ 
ing on the line, and the emission region may be extended in radius. However, 
photo-ionization calculations have shown that the line emissivities of many lines 
peak sharply in a relatively narrow range of density and ionization parameter 
(Baldwin et al. 1995). Depending on the radial variation of the density this may 
imply that the emission of a given line is localized in a relatively narrow shell 
around the continuum source. Given these results, the single radius model may 
be a reasonable approximation for a particular emission line. The physical condi¬ 
tions in the line-emitting gas are largely determined by the ionization parameter 
(the ratio of ionizing photon density to the electron density. He, e.g. Netzer 1990) 
U = Qion /where the ionizing photon flux (number of ionizing photons 
per unit time) is Qion = is fh® luminosity per unit energy, and 

Eq is the ionization potential of the line under consideration (1 Rydberg=13.6 eV 
for the hydrogen lines). Defining the ionizing luminosity. Lion = ^{^)dE 

and the mean energy of an ionizing photon, Eion = LionlQion, the BLR radius 
may be written as 


r = 


Li 


ATTcEirmUn. 


1/2 


= 13.4 


/ L44 \ 
\Unioe) 


1/2 


light — days 


( 5 ) 


where nio = ne/ 10 ^° cm“^, L 44 = Ljon/lO*^^ erg s“^, and e = Eion/Eo is the 
mean photon energy in units of £" 0 . Combining equations and |5| we get the 
expression for the black hole mass (in units of W^Mq): 


where V 3 = v{FWHM)/lO^ km s“^. Analyses of the broad emission lines in 
various AGN indicate that typical values in the gas emitting the high excitation 
broad lines are U ~ 0.1-1 and rig ~ lO^^-lO^^ cm“^ (e.g. Rees, Netzer & 
Ferland, 1989), so that Unio ~ 0.1-10. For a given ionizing spectrum, one can 
explicitly calculate the photon flux Qion, or alternatively e and Lion in eq. |5[ 
The ionizing flux is dominated by the EUV continuum in the 1-10 Rydberg 
regime, where most of the ionizing photons are emitted. Since the continuum 
in this range cannot be observed directly, it may be estimated by extrapolation 
from the nearest observable energy bands: the UV and soft X-rays. The far UV 
spectrum has been observed beyond the Lyman limit for about 100 quasars to 
wavelengths of 600 A (Zheng et al. 1996). Wandel (1997) uses an estimate of 
the ionizing flux based on the soft X-ray continuum spectrum. This approach 
is particularly appropriate for lines with high ionization potentials such as C iv 
and O VI, which are derived from more energetic photons, thus closer to the soft 
X-ray band. 


3.2. The emission-volume method 

A different approach that does not require the thin shell assumption tries to 
estimate the size of the emitting volume. Assuming this volume to be roughly 
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( 7 ) 


spherical, the line luminosity would be (Dibai 1981) 


Lline = 


where jum is the line volume emissivity and is the volume filling factor. Dibai 
(1981) took = 0.001 which is arbitrary and probably an overestimate. Wandel 
&: Yahil (1985) elaborated the method by expressing in terms of the better 
known fa, the angular covering factor: ~ where N is the column density 

of the line emitting gas. Substituting this into equation ^ gives (for the H/3 line) 

^ ^ ( lO4^^fgi-0 (^10^23^)“^/' light - days, (8) 


where N 23 = cm 2 . The implied kinematic black hole mass is obtained 

by combining this with eq. 


Mg = 0.4u| 


U2{Hf5) 

nioA^23/a/0.1 


1/2 


(9) 


Applying this method for a sample of about 90 low-redshift AGN (mainly 
Seyfert 1 galaxies and quasars), Wandel &: Yahil (1985) found a tight corre¬ 
lation between the black hole mass (inferred from the H/3 line) and the absolute 
magnitude (fig. 1). A similar value may be obtained from the sample of Joly 
et al. (1985). 



Figure 1. Optical luminosity versus the mass calculated with the emission vol¬ 
ume method from the H/3 line. 


4 






Although in part this correlation is expected, because of the known strong 
correlation between the line and continuum luminosities, the correlation found is 
actually steeper and stronger than the expected one. If true, this correlation may 
imply a universal mass-luminosity relation, which, in terms of the Eddington 
luminosity (AEdd = 1-3 x 10^^Mg erg s“^), can be expressed as logL/LEdd = 
-2 ±0.5. 


3.3. The narrow emission lines 

Similar methods can be applied to the narrow emission-line region (NLR), pro¬ 
vided the gas there is gravitationally bound. Typical parameter values for the 
NLR are UniQ = 10“^° and fa = 0.01 (e.g. Brotherton et al. 1994). Of course, 
within the distance of the NLR from the center, the stellar mass of the host 
galaxy could be comparable to the mass of the central black hole, but they 
may be correlated (see Laor 1998). Wandel Mushotzky (1986) find an ex¬ 
cellent correlation between the mass estimated from the narrow-line width and 
the black hole mass estimated from X-ray variability (fig. 3 in section 4 below). 
Wandel &: Mushotzky (1986) also apply this method to a sample of Seyfert 2 
galaxies hnding that Seyfert 2s as a group obey the same correlation, provided 
we observe only a fraction (0.1) of the emitted X-ray continuum. This is consis¬ 
tent with presence of an obscuring torus, as suggested by the unified scheme of 
AGN. 

3.4. Reverberation mapping 

Variability in the continuum UV luminosity is observed to be followed by cor¬ 
responding changes in the line luminosity with delays ranging from days in 
low-luminosity AGN to months in luminous quasars. Gross-correlating the con¬ 
tinuum and emission-line light curves gives the peak of the power in the delay 
and the characteristic size {Rr = cr) of the emission-line region. Using this size 
in eq. Q gives an alternative method to calculate the black hole mass. Gombin- 
ing all the objects with reliable reverberation data, it appears that R oc 
(Peterson 1994; Kaspi 1997). With this dependence in mind, eq. § gives a sim¬ 
ilar empirical dependence to that in eqs. ^ and M (X The effective 

BLR radii for the reverberation and ionization-parameter methods (both for the 
H/3 line) have been found to agree very well, particularly when the shape of the 
ionizing continuum is taken into account (Wandel 1997). So do of course the 
inferred masses. 

3.5. The variable component of the line profile: RMS spectrnm 

Although giving a more direct measurement of the BLR size than the photo¬ 
ionization methods, the reverberation mapping kinematic method has two draw¬ 
backs: 

• Very few AGN (about a dozen) have reliable reverberation data. Mea¬ 
suring the reverberation radius requires many observations and is a large 
project. 

• It is not clear that the line variability size measured by the technique of 
reverberation (or echo) mapping is the same one contributing the emission¬ 
line profile. 
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The line profile is given by 

P Rout 

/(AA) =4tt 6[v{R) sm(i) - c{AX/X)]E{R)R^dR 

Rin 

where E{R) is the volume line emissivity at a radius R and i is the inclination 
angle to the line of sight. This may be used to define a velocity-weighted radius, 
e.g. Rfwhm, associated with the FWHM of the line profile. A related radius 
is the emissivity-weighted radius, that may be defined by 

pRout pRout 

Rem = / E{R)R^dR / / E{R)R^dR. (11) 

The reverberation radius is obtained by cross correlating the line and continuum 
light curves. 

For a thick geometry these radii may be quite different. One solution to 
this difference is to combine the reverberation-mapping radius with the velocity 
inferred from the variable component of the line profile. In this method, the 
distance of the line-emitting gas from the cetral source (the reververation radius) 
and its velocity dispersion (given by the the rms spectrum, which reflects only the 
variable part of the emission-line profile) both refer to the same gas. Comparing 
the rms reverberation and the ionization-parameter methods, one may calibrate 
the latter, which is easier to use and may be applied also for objects without 
variability data. 

3.6. Uncertainties and potential error sonrces 

There are several drawbacks which introduce possible errors into the emission¬ 
line kinematic methods: 

• Uncertainty and scatter in the parameters. In the ionization method we 
do not know the values of n and U. 

• In the volume-emission method the parameters N, fa and n are uncertain. 

• The velocity of the emitting gas may be anisotropic. 

• The illumination of the clouds may be anisotropic. 

• The motions may be due to non-gravitational forces. 


4. Time variability 

The AGN continuum is known to vary on time scales that depend on the wave¬ 
length and luminosity. While for the UV the time scale ranges from days to 
months (depending on the luminosity), the time scale for variations in X-rays 
is of the order of hours. As we have seen above, the variability can be used to 
place a lower limit on the black hole mass. If we knew the radial dependence of 
the local X-ray spectrum in the inner disk, the effective emission radius in eq. 
would be 

pRout pRout 

Reff = / Fa,{R)R^dR / / E^{R)RdR (12) 

Rin Rin 
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where Fx{R) is the X-ray emissivity per unit area from the disk at radius R. 
Since there is no widely accepted model for the X-ray continuum emission in 
AGN, we assume Reff is a few Schwarzshild radii, e.g. where the emission 
from the thin disk around a Schwarzshild black hole peaks. For this value eq. |l| 
gives Mg < 0.5At/10^ s. Barr & Mushotzky (1986; hg. 2) have found a good 
correlation between At and the X-ray luminosity. 



log L(x) 2-10 kev (ergs/sec) 


Figure 2. X-ray variability time scale versus X-ray luminosity. Circles — Barr 
&: Mushotzky (1986), squares — EXOSAT data, crosses — radio-bright objects. 
Solid lines are best fits (with and without Galactic black holes). The dashed 
line indicates the Eddington limit assuming R = 5Rs. 

While the variability masses of quasars and Seyfert Is seem to give Edding¬ 
ton ratios well below unity, some bright radio objects (e.g. BL Lacs) as well 
as some Narrow-Line Seyfert 1 galaxies (NLSl) have values close to unity and 
even larger than unity, e.g. PHL 1092 (see Brandt &: Boiler 1998). Taking into 
account the data on Galactic black holes and bright radio sources, one gets a 
smaller slope in the L-At plane, which may be interpreted (assuming M oc At) 
as the L/M ratio increasing with luminosity. It may therefore be of interest to 
compare the masses derived by the variability method (which admittedly gives 
only an upper limit) with the predictions of other methods, in order to find how 
tight the upper limits are, or how well they correlate with the real mass. Wan- 
del & Mushotzky (1986) find a very good linear correlation between the X-ray 
variability time and the kinematic mass (calculated from the [O iii] line; see fig. 
3), which supports the hypothesis that At (X-ray) and the FWHM of [O iii] 
both are closely related to the central mass. Interestingly, in the sample used by 
Wandel & Mushotzky (1986), the correlation of At with the mass derived from 
H/3 is less good than that with the mass derived from [O iii]. For broad lines, 
a good correlation between At (X-ray) and the H/3 line width is obtained in a 
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mixed sample containing rapidly variable NLSl and normal Seyfert Is (Wandel 
& Boiler 1998). 


log M(At) 

678 9 10 



Figure 3. X-ray variability time versus [O iii] kinematic mass estimate (Wandel 
and Mushotzky 1986). 

5. Accretion-disk modeling 

Many authors have tried to fit accretion disk spectra to the observed AGN 
continuum, thus finding the accretion disk parameters (M, M, a and the angular 
momentum of the black hole) which best fit the observed continuum in the UV 
(e.g. Wandel & Petrosian 1988; Sun & Malkan 1991) or in the soft X-rays (e.g. 
Laor 1990). 

5.1. The accretion disk - black body spectral relation 

The thin accretion disk spectrum is a multi-temperature black body one, and 
trying to fit the UV bump (or the soft X-ray spectrum) is essentially equivalent 
to using eq. ^ in order to find the temperature-mass: Mg « /T^{r/\{)) 











where r = R/Rg. The second accretion disk parameter, M, is determined by 
the observed luminosity, schematically via the relation L = eM(P‘. An alter¬ 
native approach relates the peak in the accretion disk spectrum to the thin 
accretion disk parameters by assuming the spectrum is dominated by black 
body emission near the maximum-surface emissivity radius. This gives 

Emax ~ 3kTbb{5Rs) = 10(m/M8)^/^ eV (13) 

where m = M/MEdd is the Eddington ratio. More sophisticated approaches, 
like taking as Tmax the temperature at the inner boundary between the black 
body and the radiation-dominated disk or the temperature at which the inner 
disk becomes optically thin give Emax ~ eV. While in the UV band 

the thin accretion disk multiple black body spectrum may be a good approxi¬ 
mation, the soft X-rays may be produced by a hotter medium due to processes 
other than black body emission, such as Comptonization, a two-temperature 
disk (Wandel & Liang 1991) or hot corona (Haardt & Maraschi 1993; Czerny, 
Witt & Zycki 1996). 



Figure 4. Accretion-disk black hole evolution tracks in the a — L plane (Wandel 
& Petrosian 1988). Crosses: Seyfert galaxies, triangles, open circles and filled 
circles: low, medium and high redshift quasars, respectively. Continuous lines: 
constant mass; dashed lines: constant M/M (note that, unlike in the text, the 
labeling of the curves in the figure follows the notation m = 16.7M/MEdd)- 

5.2. Deriving M and M from the UV spectrum 

It is possible to perform a more refined treatment, calculating for each pair of 
accretion disk parameters M and M not the total luminosity and blue-bump 
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temperature, but the actual observables, e.g. the UV luminosity and spectral 
index. Wandel & Petrosian (1988) have calculated the accretion disk flux and 
spectral index at 1450 A for a grid of accretion disk parameters (M, M). In¬ 
verting the grid they have obtained contours of constant black hole mass and 
constant Eddington ratio (rh) in the L — ajjv plane (fig 4). Plotting in this 
plane samples of AGN it is possible to read from the contours the correspond¬ 
ing accretion disk parameters. Comparing several groups of AGN a systematic 
trend appears: higher redshift and more luminous objects tend to have larger 
black hole masses and luminosities closer to the Eddington limit (see table 1). 
Similar results are obtained by Sun & Malkan (1991). 


Table 1. Grouping of AGN accretion disk parameters. 

AGN group Log F1450 Log M/M© M/MEdd 


Seyfert galaxies 

Low z quasars 
Medium 2; quasars 
High 2; quasars 

28- 29.5 

29- 30.5 

30- 31.5 
31-32 

7.5-8.5 0.01-0.5 

8-9 0.02-0.1 

8- 9.5 0.1-0.5 

9- 9.5 0.03-2 

Table 2. L/M ratios found by various 

methods. 

Method 

L/ L-^dd 

Reference 

Kinematic 

Emission Volume (H/3) 
Ionization parameter (H/3) 

IP -|- Reverberation(H/3) 

0.003-0.03 

0.01-0.1 

0.01-1 

Wandel & Yahil (1995) 

Joly et al. (1995) 

Wandel (1997) 

X-ray variability 

X-rays 

X-rays-l-OIII photoionization 
Soft X-rays, NLSl 

0.003-0.1 

0.003-0.1 

0.001-0.1 

Barr & Mushotzky (1986) 
Wandel & Mushotzky (1986) 
Wandel & Boiler (1998) 

Accretion-disk fitting 

Modified BB disk 

Kerr black hole disk 

GR effects 

Disk-I-Corona 

0.001-2 

0.01-2 

0.003-0.3 

0.001-0.3 

Wandel & Petrosian (1986) 
Sun & Malkan (1991) 

Laor (1990), Laor (1993) 
Czerny, Witt & Zycki (1996) 
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6. Summary 


The L/M ratios found by various methods are summarized in table 2. All 
emission-line kinematic methods give a relatively narrow L/M range over a large 
range in luminosity, suggesting a universal L/M relation in AGN. However, this 
may be related to the fact that the samples with available data are mostly low 
redshift, relatively low-luminosity objects. The two other methods — X-ray 
variability and accretion disk fitting of the UV bump — suggest a trend of the 
Eddington ratio increasing with luminosity. 
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Discussion 

Brad Peterson: I believe that a fairly reliable virial mass can be obtained from re¬ 
verberation mapping data alone. The continuum emission-line cross-correlation 
function gives a centroid that measures the responsivity-weighted radius of the 
BLR. The spectra can be combined to form mean and rms spectra. The rms 
spectrum identifies the variable part of the emission line, and the line width 
used for virial estimates should be that measured from the rms spectrum. This 
way both the characteristic line width and scale size refer to the same gas, that 
which is varying at the time of the experiment. 

Amri Wandel: I fully agree. This method actually gives “the best of both 
worlds”; a direct measurement of the distance along with the appropriately 
weighted velocity. It is interesting to compare the rms reverberation and the 
ionization-parameter methods, in order to calibrate the latter, which is used 
more readily and also for objects without variability data. 

Miehael Corbin: I wish to comment that there exists a strong correlation be¬ 
tween the asymmetry of the G iv A1549 and the H/3 profiles and the luminosity 
of the soft X-ray and near UV continuum, such that as luminosity increases, the 
profiles become more red asymmetric. I have been able to model these redward 
asymmetries as the effect of gravitational redshift, which is proportional to the 
mass of the black hole. This correlation is therefore consistent with a universal 
black hole mass-luminosity relationship. 

Martin Gaskell: It is interesting that the accretion disk method shows M tending 
to MEdd when we go to high luminosities. From reverberation mapping (Ko- 
ratkar & Gaskell 1991, ApJL) we found 3G273 near MEdd- Prab Gondhalekar 
has found this to be generally true for high luminosity quasars. However, from 
our analyses of 3G273, I think that beaming of the continuum is causing an 
underestimate of the radius in reverberation mapping. 

Amri Wandel: The trend of the L/M ratio to be larger for more luminous objects 
appears to be associated only with the accretion disk fitting method of mass 
estimation. When the mass is estimated by the ionization-parameter method, 
this trend is much weaker or even consistent with M oc L (hence L/M ~ const). 
The correlation between L/M and the luminosity in the accretion disk method 
may therefore reflect a bias related to the thin accretion disk equations or the 
UV-X-ray continuum of AGN. 

Jean Clavel: Observationally, how does one define a characteristic “fastest” 
X-ray variability? The power-density has a slope close to unity (e.g. is self¬ 
invariant), so that there is no characteristic time scale and the probability of 
observing an event of AL within a time At increases with the length of the 
observation. 
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Amri Wandel: The variability time scale has been determined only from large 
events (change of the flux by a factor of 1.5 or more), and the characteristic time 
has then been taken as the e-folding time of that change. 
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